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Abstract The parabolic Anderson problem is the Cauchy problem for the heat 
equation d t u(t,z) = Au(t ,z) + £(t ,z)u(t ,z) on (0,°°) x 1r with random potential 
(t,z) ■ z€ Z d ) and localized initial condition. In this paper we consider potentials 
which are constant in time and independent exponentially distributed in space. We 
study the growth rate of the total mass of the solution in terms of weak and almost 
sure limit theorems, and the spatial spread of the mass in terms of a scaling limit 
theorem. The latter result shows that in this case, just like in the case of heavy tailed 
potentials, the mass gets trapped in a single relevant island with high probability. 

1 Introduction and main results 
1.1 Overview and background 

We consider the heat equation with random potential on the integer lattice Z d and 
study the Cauchy problem with localised initial datum, 

d,u(t,z) = Au(t,z) + £,(t,z)u(t,z), for (t,z) G (0,°°) x If' , 

lim«(f,z) = l (z), for zeZ d , 

40 

where 

(4/)(z)=£[/W-/(z)], forzeZ f/ ,/: 1 d ->R, 

is the discrete Laplacian, and the potential (t,z) : t > 0,z <G I d ) is a random field. 
This equation is known as the parabolic Anderson model. 
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In the present paper we assume that the potential field is constant in time and 
independent, identically distributed in space according to some nondegenerate dis- 
tribution. Under this hypothesis the solutions are believed to exhibit intermittency, 
which roughly speaking means that at any late time the solution is concentrated in a 
small number of relevant islands at large distance from each other, such that the di- 
ameter of each island is much smaller than this distance, see Figure 1 for a schematic 
picture. The relevant islands are located in areas where the potential has favourable 
properties, e.g. a high density of large potential values. As time progresses new rel- 
evant islands emerge in locations further and further away from the origin at places 
where the potential is more and more favourable, while old islands lose their rele- 
vance. The main aim of the extensive research in this model, which was initiated by 
Gartner and Molchanov in [3, 4], is to get a better understanding of the phenomenon 
of intermittency for various choices of potentials. 




Fig. 1 A schematic picture of intermittency: the mass of the solution is concentrated on relevant 
islands (indicated by shaded balls) with radius of order a, and distances of order r, 3> a t . 

Natural questions about the nature of intermittency are the following: 

• What is the diameter of the relevant islands? Are they growing in time? 

• How much mass is concentrated in a relevant island? 
How big is the potential on a relevant island? 

• Where are the relevant islands located? What is the distance of different islands? 

• How many relevant islands are there? 

• How do new relevant islands emerge? What is the lifetime of a relevant island? 
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Explicit answers to these questions and, more generally, results on the precise 
geometry of solutions to the parabolic Anderson model are typically very difficult 
to obtain. In the related context of Brownian motion among Poissonian obstacles, 
Sznitman [14] provides methodology to study properties of Brownian paths con- 
ditioned on survival, which offer a possible route to the geometry of solutions, 
at least in the case of bounded potentials. In a seminal paper Gartner, Konig and 
Molchanov [2] follow a different route to analyse size and position of relevant is- 
lands in the case of double exponential potentials. Their results also offer some 
insight into potentials with heavier tails. In [7] and [9] a complete picture of the 
geometry of the solutions is given in the case of Pareto distributed potentials, build- 
ing on the work of [2]. In this case of an extremely heavy tailed potential it can 
be shown that, for any e > at sufficiently late times, there exists a single point 
carrying a proportion of mass exceeding 1 — £ with probability converging to one. 
This point constitutes the single relevant island and very precise results about the 
location, lifetime and dynamics of this island can be obtained, see also [10] for a 
survey of this research. 

For more complicated potentials however, one has to rely on less explicit results. 
A natural way forward is to investigate the growth rates of the total mass 

U(t):= £ u(t,z) 

zez d 

of the solution. If the potential is bounded from above we define the (quenched) 

Lyapunov exponent as 

A := limL, where L, := -logt/(f), 

whenever this limit exists in the almost sure sense. If the potential is unbounded one 
expects superexponential growth and is interested in an asymptotic expansion of L t . 
If the tails of the potential distribution are sufficiently light so that the logarithmic 
moment generating function 

H(x) :=log£e^ (0) 

is finite for all x > 0, a large deviation heuristics suggests that, we get 

H(B ar d ) 1 
L, = j(k + o(1)), almost surely as 1 1 °°, 

where a,j3 are deterministic scale functions and K is a deterministic constant. Ac- 
cording to the heuristics, the quantity a, can be interpreted as the diameter of the 
relevant islands at time t, and the leading term as the size of the potential values on 
the island. The constant k is given in terms of a variational problem whose max- 
imiser describes the shape of a vertically shifted and rescaled potential on an island. 
More details and a classification of light-tailed potentials according to this paradigm 
are given in [5]. 
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If the potential is such that the moment generating functions do not always exist, 
this approach breaks down. Indeed, one can no longer expect the leading terms in 
an expansion of L t to be deterministic. Instead, one should expect the solutions 
to be concentrated in islands consisting of single sites and the expansion of L t to 
reflect fluctuations in the size of the potential on these sites. One would expect the 
sites of the islands to be those with the largest potential in some time-dependent 
centred box and the fluctuations to be similar to those seen in the order statistics 
of independent random variables. This programme is carried out in detail in [6] for 
potentials with Weibull (stretched exponential) and Pareto (polynomial) tails. In the 
present paper we add the case of standard exponential potentials and present weak 
(see Theorem 1) and almost sure (see Theorem 2) asymptotic expansions for Lt in 
this case. These results are taken from the first author's unpublished master thesis [8] 
and were announced without proof in [6]. 

Very little has been done so far to get a precise understanding of the number and 
position of the relevant islands, the very fine results for the Pareto case being the 
only exception. A natural idea to approach this with somewhat softer techniques is 
to prove a scaling limit theorem. To this end we define a probability distribution v t 
on 1r associating to each site z a weight proportional to the solution u(t,z), i.e. 

v ' : = E ^S(z),foranyf>0, 

zez d U{ ) 

where 8(z) denotes the Dirac measure concentrated at z € M. d . For a > 0, we also 
define the distribution of mass at the time t in the scale a as 

zez' 1 u rJ 

which is considered as an element of the space ^(R d ) of probability measures 
on W 1 . Identifying the scale r t of the distances between the islands and the origin, 
intermittency would imply that islands are contracted to points and that vp con- 
verges in law to a random probability measure, which is purely atomic with atoms 
representing intermittent islands and their weights representing the proportion of 
mass on the islands. In the case of Pareto potentials such a result follows easily 
from the detailed geometric picture, see [9, Proposition 1.4], but in principle could 
be obtained from softer arguments. It therefore seems viable that scaling limit the- 
orems like the above can be obtained for a large class of potentials including some 
which are harder to analyse because they have much lighter tails. 

In Theorem 3 of the present paper we show that in the case of exponential poten- 
tials for r t =t/ log log t the random probability measures v/' converge in distribution 
to a point mass in a nonzero random point. In particular this shows that for expo- 
nential potential we also have only one relevant island. Moreover, the solution of 
the parabolic Anderson problem spreads sublinearly in space. Our arguments can 
be adapted to the easier case of Weibull, or stretched exponential, potentials, where 
there is also only one relevant island but the solution has a superlinear spread. These 
results are new and open up possibilities for further research projects, which we 
briefly mention in our concluding remarks. 
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1.2 Statement of results 

We now assume that z <E Z d ) is a family of independent random variables 

with 

P(§(z) >x) = e~*forx>0. 

Suppose (u(t,z) : t > 0,z <G Z ) is the unique nonnegative solution to the parabolic 
Anderson model with this potential, and let (17(0 • t > 0) be the total mass of the 
solution. We recall that 

L t = jlogU(t) 

and first ask for a weak expansion of L t up to the first nondegenerate random term. 
This turns out to be the third term in the expansion, which is of constant order. In 
the following we use =>■ to indicate convergence in distribution. 

Theorem 1 (Weak asymptotics for the growth rate of the total mass). 

We have 

L, — a! log f + if log log log f =>X, 
where X has a Gumbel distribution 

P(X <x) =exp{ -2 d e~ x+2d } forxtR. 

In an almost sure expansion already the second term exhibits fluctuations. 

Theorem 2 (Almost sure asymptotics for the growth rate of the total mass). 

Almost surely, 

Lt-dlogt 

hm sup = 1 , 

ff. log log t 

and 

. Lt-dlogt 

hminf- — - — = -(d+1). 

ft~ log log log t 

Remark 1. Note that neither of these almost sure asymptotics agree with the asymp- 
totics 

Lt-dlogt . 

hm = —d in probability, 

<t~ log log log t 

which follows from Theorem 1 . The almost sure results pick up fluctuations on both 
sides of the second term in the weak expansion, with those above being signifi- 
cantly stronger than those below the mean. This is different in the stretched expo- 
nential case studied in [6], where the liminf behaviour coincides with the weak limit 
behaviour. The limsup behaviour in the exponential case is included in the results 
of [6] and therefore not proved here. 
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Recall that the distribution of the mass of the solution at time t > and on the 
scale a > is defined as a (random) element of the space ^(W 1 ) of probability 
measures on M by 

zEZ d [ > 

The following theorem is the main result of this paper. 

Theorem 3 (Scaling limit theorem). Defining the sublinear scale function 

t 

r ' log log f' 

we have 

limv/ 7 = 8{Y) in distribution, 

where 8 (x) denotes the Dirac measure concentrated in x £ Mr 1 and Y is a random 
variable in M with independent coordinates given by standard exponential vari- 
ables with uniform random sign. 

Remark 2. In the case of a Weibull potential with parameter < y < 1 given by 

P{%{z) >x) =e~ x7 forx>0, 
a variant of the proof gives convergence of vp for the superballistic scale function 

log log t 

to a limit measure 8(Y) where the components of Y are independent exponentially 
distributed with parameter d 1 ^ 1 ^ and uniform sign. Details are left to the reader. 



2 Proof of the main results 



2.1 Overview 



The proofs are based on the Feynman-Kac formula 

n(f,z)=E exp{J^(X s )ds}l{X, = z} 

where (X s : s > 0) is a continuous-time simple random walk on 7L d started at the 
origin and the probability P and expectation E refer only to this walk and not to the 
potentials. Recall that (X s : s > 0) is the Markov process generated by the discrete 
Laplacian A featuring in the parabolic Anderson problem. It is shown in [3] that the 
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Feynman-Kac formula gives the unique solution to the parabolic Anderson problem 
under a moment condition on the potential, which is satisfied in the exponential case. 
By summing over all sites the Feynman-Kac formula implies that the total mass is 
given by 

C/(0=E[exp{y %(X s )ds} . 

An analysis of this formula allows us to approximate L t = |logt/(?) almost surely 
from above and below by variational problems for the potential. These variational 
problems have the structure that one optimizes over all sites z g Z d the difference 
between the potential value t, (z), corresponding to the reward for spending time in 
the site, and a term corresponding to the cost of getting to the site, which is going to 
infinity when z — > 00 and thus ensure that the problem is well-defined. 

We can use the result for the lower bound given in [6, Lemmas 2.1 and 2.3]. Here 
and throughout this paper we use | ■ | to denote the ^'-norm on IR''. 

Lemma 1 (Lower bound on L t ). Let 

Mf):=max{^(z)-Mlog^(z)}, 

then, almost surely, for all sufficiently large t, we have 

Lt >N(t)-2d + o(l). 

The appearance of % (z) in the cost term can be explained by the fact that part 
of the cost arises from the fact that the optimal paths leading to z spend a positive 
proportion of the overall time traveling to the site and therefore miss out on the 
optimal potential value for some considerable time, see Section 1.3 in [6] for a 
heuristic derivation of this formula. 

The corresponding upper bound will be our main concern here. 
Lemma 2 (Upper bound on L,). For any c > let 

N c (t):= max {§(*) - M (loglog|z| +c)\. 

f/(logf) 2 <]z|<rlogf I. f > 

Then, for any £ > there exists c = c(s) > such that, almost surely, for all suffi- 
ciently large t, we have 

Lt <N c ( e) (t)-2d + e+o(l). 

This lemma will be proved in two steps: We first remove paths that do not make 
an essential contribution from the average in the Feynman-Kac formula using an 
ad-hoc approach, see Lemma 7 and Lemma 8. Then we use the properties of the re- 
maining paths to refine the argument and get an improved bound, see Proposition 1 . 
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The variational problems for the upper and lower bound can then be studied using 
an extreme value analysis, which follows along the lines of [6]. It turns out that the 
weak and almost sure asymptotics of the two problems coincide up to the accuracy 
required to prove Theorem 1 and Theorem 2. 

For the proof of the scaling limit we need to give an upper bound on the growth 
rate of the contribution of all those paths ending in a site at distance more than 8r t , 
for some 8 > 0, from the site with the largest potential among those sites that can 
be reached by some path with the same number of jumps. This bound needs to be 
strictly better than the lower bound on the overall growth rate. To this end, in a first 
step, we again use Lemma 7 and Lemma 8 to eliminate some paths using ad-hoc 
arguments. In the second step we remove paths that never hit the site with largest 
potential that is within their reach. This is done on the basis of the gap between 
the largest and the second largest value for the variational problem in the upper 
bound. In the third step it remains to analyse the contribution of paths that hit the 
optimal site but then move away by more than 8r,. Again it turns out that the rate 
of growth of the contribution of these paths is strictly smaller than the lower bound 
on the growth rate of the total mass. Proposition 1 is set up in such a way that it can 
deal with both the second and third step. We conclude from this that the solution is 
concentrated in a single island of diameter at most 8r t around the optimal site. An 
extreme value analysis characterizes the location of the optimal site and concludes 
the proof of Theorem 3. 

The remainder of the paper is structured as follows: In Section 2.2 we give some 
notation and collect auxiliary results from [6]. Section 2.3 contains the required up- 
per bounds and constitutes the core of the proof. Section 2.4 studies the variational 
problem arising in the upper bound. Using these approximations we complete the 
proof of Theorem 2 in Section 2.5 and of Theorem 1 in Section 2.6. The proof of 
the scaling limit theorem, Theorem 3, is completed in Section 2.7. 



2.2 Auxiliary results 

Let By = {|z| < r} be the ball of radius r centered at the origin in Z . The number l r 
of points in B r grows asymptotically like r d . More precisely, there exists a constant 
K d such that, lim, .^oo l r r~ d = K d . We define M r = maxui< r ^(z) to be the maximal 
value of the potential on B r . The behavior of M r is described quite accurately in [6, 
Lemma 4.1], which we restate now. 

Lemma 3 (Bounds for M r ). Let 8 € (0, 1) and c> 0. Then, almost surely, 

M r < Jlogr + loglogr+ (log log r) 5 for all sufficiently large r, 
M r > dlogr — (1 + c) logloglog r for all sufficiently large r. 

In particular, for any pair of constant c\ and cj_ satisfying c\ < d < cj, we have 

c \ l°g r < M r < C2 log r for all sufficiently large r. 
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Let My ] denote the z'-th biggest value taken by the potential in the ball of radius r 
centered at the origin. The next lemma gives us estimates for upper order statistics 
for the potential. 

Lemma 4 (Rough asymptotic behaviour for upper order statistics). Let < /3 < 

1 be a fixed constant. Then, almost surely, 



lim 



M„ 



(L«0J) 



log« 



= d-p. 



Proof. Recalling that /„ is the number of points in a ball of radius n in Z d we get 

L 

i=0 



L«^J - i /■ \ 

•(Mf^<jcj= £ ^y- x '(\-e- x )'"-' . (1) 



We fix £ > and infer that 



< (n^ + l) ((K cl + o(l))nP +£ y exp[-(^ + o(l))n' ,+£ 

-nP + \K d + 0{\)) 



exp 



Since this sequence is summable, we can use the Borel-Cantelli lemma to obtain 
the lower bound. Similarly, for the upper bound, we use (1) to get 



p(M ( y^ ] > (</-/J + e)logn) < £ 



,-(d-p+e)i 



(2) 



We now use a rough approximation for the binomial coefficient, namely 



when i is big enough. Combining this with (2) and using that the first term in the 
ensuing sum is the largest, we obtain, for all sufficiently large n, 



In 



p(M^^ > {d-p + e)logn) < £ 



el„ 



<ln 



'In 



■.d+e 



< 



Using the Borel-Cantelli lemma again we obtain an upper bound, completing the 
proof of our statement. 
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Let 0<(7<p<jbe some fixed constants. We define 

k n = \n a J and m„ = \n p \ 

Combining Lemma 3 and Lemma 4, we get the following result. 
Lemma 5. For any constant c > 0,for all sufficiently large n, we have 

(i) M^-M^ > (a-c)logn; 

(ii) M„" ] - M ( n" ] >{p-a-c) logn. 

Finally, we use Lemma 3 to give a lower bound for N_(t). 
Lemma 6 (Eventual lower bound for N(t)). For any small e > 0, we have 

N(t) > c/logf- (d + 1 +e)logloglogf, 

for all sufficiently large t, almost surely. 

Proof. Using Lemma 3 we get, for any fixed c > and C2 > d, 

r r r 

N_(t)>max dlogr — (1 +c)logloglogr loglogr logC2 , 

r>o L t t i 

if the maximum of the expression in the square brackets (which we denote by f(r)) 
is attained at a point r t , large enough so that Lemma 3 holds. 
The solution r = r, of f\ (r) = satisfies 

d = loglogr 
r t 

Writing r, = t(p(r t ), where (p(r) — c/(loglogr) _1 (l +o(l)) we get that 

log<p(r) = -logloglogr + logc/ + o(l) (3) 

and hence logr, = logf + log(p(r t ) = logf + o(logr r ), which implies logf>/logf = 
1 +o(l). Note that this implies r t — > °° as t — > «, which justifies a posteriori the 
application of Lemma 3. Combining this with (3) we get, 

f(r t ) =d(log(tq>(r t ))) - (1 +c) log log logr, - 9 (r r ) (log logr, +logc 2 ) 
= <ilogr- (1 +<i + c) log log logf + 0(1). 



2.3 Upper bounds 



We start by showing ad-hoc bounds for the growth rates of the contribution of certain 
families of paths. These can be compared to the lower bound for the growth rate of 
U(t) showing that the paths can be be neglected. For a path (X s : s > 0) on the 
lattice Z d we denote by J, the number number of jumps up to time t . Recall that 
M« denotes the k th largest potential value on the sites z G If 1 with |z| < «. 
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Lemma 7. Fix < a < A and k n = n a . Let 

tf 2 (t)=E[exp{ jT§&)dj}l{^ <J t <flogf,max^(X s ) << J ''} 
77ze« 

hm - log — — = — °°. 
rt°° t B U{t) 

Proof. Simply replacing | in the integral by the maximum, we get 

U 2 (t)= Y E[exp( / t(X s )ds)t\j t =n,max£(X s ) <M { * n) } 
t/(lo g t)Hn<tlo g t 1 Uo J 

< £ e ^"'p(y,=«)< max X'"'. 
(/(logr) 2 <»<rlogf tl(\ogtf-<n<t\ 0% t 

By Lemma 4 we have M„ = (d — (?) logn + o(logn) and hence 

j\ogU 2 (t) < (d-o)logt + o(\ogt), 

so that the result follows by comparison with Lemma 1 and Lemma 6. 
Lemma 8. Let 

exp { \\ (X s ) ds] (l{J, < ^ } + l{J t > dog/} 



U 3 (t)=E 
Then 

y 1 1/3(0 

hm - log — — = — oo. 
Proof. We first show that almost surely, 

ilogC/ 3 (0< max {M B -^log-^-} -2J + o(l). 

Indeed, we have 

(2dt) n e- 2 '" 



U 3 (t)< E e ' M "P(/,=n)= £ e 

{*<f/(logr) 2 } {„<r/(logr) 2 } 
U{n>flogi} U{nxlog(} 



tM„ 



To estimate n\ we use Stirling's formula, 

n\ = V27m(-)" e S(n \ with limS(n) =0. 

V e / nt°o 



(4) 



< V exy(tM n -2dt + nlog2dt -lognl) . (5) 

{<.<r/(logf) 2 } 
U{/j>rlogr} 
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Fixing some e > we know from Lemma 3, that M n < (d + e)logn for all suffi- 
ciently large n, so for f large enough, we obtain for all n > f logf, 

tM„ - 2dt + n log 2dt - log n ! < t (d + e ) log n - n log ^ - 8 (n) 

<^ + e )log„(l-ig)log(^) +0 (l)) 
< — 21og«, 

by noticing that n i-> log is decreasing on (t logf ,«>). Hence, almost surely, 
exp(tM n -2dt + nlog2dt-lognl) =o(l), 

»>rlogr 

so that using (5) the following upper bound for t/3 

U3U) < — 1 —~ max exp(/M„-2df + nlog2d/-logrt!) + o(l) 
(logf) 2 «<f/(logr)2 



< 



f 



max exp (tM n - 2dt - n log ^ + o (/)) + o( 1 ) 



(logf) 2 «<f7(iog*) 2 
and hence (4) follows. As a second step we show that 

-logE/ 3 (*) < dlogt- (2d- l)loglogf + o(loglogf). (6) 

Recall that r H> M r is a non-decreasing function and check that 
r r 

ri — > -log^i — is decreasing on (0,2det), 
t 2det 

hence, replacing r in the bracket by f/(logf ) 2 

1 = M r/(log0 2 +°( 1 )- 



max 

r<t/(\ogt) 2 



M '-7 log 2dVt 



By Lemma 3 we have M r < d log r + log log r + o (log log r) for all sufficiently large r, 
we get, for f large enough 



max 

r<r/(logr) 2 



M r --\og-^— < dlogt- (2d- 1) log logf + o(loglogf), (7) 



f 2c/ef 

and combining (4) and (7), we have proved (6). Using Lemma 1 and Lemma 6, 

ylog^|| < l -logUi(t)-N(t)-2d + o(l) 

< — (2d— l)loglogf + o(loglogf) — » —00, 
and hence our statement is proved. 
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The following versatile upper bound is the main tool in the proof of all our the- 
orems and will be used repeatedly. Note for example that, together with Lemmas 7 
and 8 it implies Lemma 2 if the parameters in (;/) are chosen as k = 1 and 8=0. 

Proposition 1. For a path (X s : s > 0) on the lattice 1 d we denote by J t the number 
of jumps up to time t. We denote by M„ the k th largest potential value on the sites 
zeZ d with \z\ < n, and let Z„ be the site where this maximum is attained. Further 
fix < tT < j and let k n = \ n° J and a t \. 0. 

(a) Forn G N let 

C/, ( " ) (f)=E[exp( f ${X s )ds\l{J, =n)l{ max§(X,) >M< W ) . 

L I Jo > 0<s<t i 

Then, for all £ > there exists C £ > such that uniformly for allta t <n <t log t, 
- log (f ) < M' 1 ' - |- ( log log n - C e ) + £ - 2d + o ( 1 ) as 1 1 °°. 

(b) For fixed 8 > and k,n£N let 

Uf' k ' n) {t)=v\^[j t ^{X s )ds) 

l{J t = n}l{ sup ^(X s )^{Mil\...,M l t i] }} 

0<s<< > 

1{Z« e {X s : < s < t }, \X t - Zf I > fir,}] . 

Then, almost surely, 



uniformly in k < k„ and t y < n < t Or, 



we have that 



i 7 M| 

■ log Uf' k ' n) (t ) < M« - i=2-L log log |Z« | - 2d - 5 + o(l) as / 1 ~. 



The first step in the proof is to integrate out the waiting times of the continuous 
time random walk paths. The following fact taken from [6] helps with this. 

Lemma 9. Let t]Q,...,T] n be fixed real numbers attaining their maximum only once, 
i.e. there is an index <k <n with rj^ > r\ifor all i ^ k. Then, for all t > 0, 

f f" _1 / " _1 \ ~i r" _1 ~i l 

/ exp<^ y>7?i+ [t-Tti)ri n \l\ Yti<t\dto...dt n -i<e tnk \\ . 

Ju% I to V to > S { to S t/kVk-ri, 
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Proof. First, we prove the result for the case k — n, i.e. rj„ > r\i for all ; < n. We 
have 

n— 1 , n— 1 ^ N , n— I 

-1 



n exp| £ '/ty + (t - £ f <) 1" } 1 { L f >' < f ' f '' ^ 0V/ < « - 1 } ^0 • • ■ 

1=0 i'=0 (=0 

= <? ,7? " / exp{ Y,ti{i]i -7] n )}t\ Yti<t}dt ...dt n -i 

/■ r" _1 ~i l 
<e">" / exp V f; (tj,-- rj n )\dt ...dt n - 1 = e"l"n ■ 

Now we show that any permutation of the indices does not change the value of the 
integral above and this will be sufficient to prove the statement. First, it is obvious 
that transposition of i and j does not change the integral if i,j < n — 1. Now we 
consider the case of a transposition x on j and n, where j < n. We change variables 
such that f- = ti if / 7^ 7, i < n — 1 and t'-=t — £"~ U, and get 

. n— 1 . n— 1 . - 



/ exp { £ f/77; + (t - £ /,) tj„)i{ £ t, < t] dt . . . a„_i 

= / exp { £ f ; '77 T(l) - (t - £ f A i} t(B) } 1 

n-l 

x j £ ?■ < f ,t'j > 0V7 < n- 1 \dt' Q . . . dt' n _ u 



which completes the proof. 

For the proof of Proposition 1 (b) denote by 

^ ( " J,) = {.v = Oo, • • • ,Vn) -yo = 0, [yt-i-yil = l, 

{yo, ■ ■ ■ ,y n ] n {Z<», . . . ,zr"} = 0,Z« G {yo, • • • ,?„}, b„ -Zf\ > 8r, } 

the set of all 'good' paths and let (t,) be a sequence of independent, exponentially 
distributed random variables with parameter 2d. 

Denote by E the expectation with respect to (t,). We have 

U['**\t) = £ (2d)-"E [exp { £ x£ (y,) + (t- £ x,) % (y„)} 

v( = ca>(8,k.n) i=Q i'=0 

yG '^ (8) 

-1 



l{£^<f,£T,->r} 



x 

1=0 i=0 
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In the further proof we apply Lemma 9 to the values of the potential ^ along a 
path y. However, to do so we need the maximum of E, along the path y to be attained 
only once. Therefore we have to modify the potential along the path slightly. 

We fix y G ,^(" J '> and let 

/(y)=min{/e{0,...,«}:y,=4 i »} 

be the index of the first instant where the maximum of the potential over the path is 
attained. Now we define a slight variation of ^ on y in the following way. Fix e > 
and define & : {0, ...,«} -> R by % = | (y,) if i ± i(y), and % {y) = £ (y ;(y) ) + e. We 

obtain, using t, (y,) < that 

. n— 1 , n— 1 v » , n— 1 



exp { £ (y ( ) + (f - £ T,) 4 (y„) jl{ £ Tj < r , £ t, > t } 

k i=0 V i=0 ' ' V 1=0 1=0 

< E [exp { £ +(t- £ t ; ) §>}l{ £ T, < t, £ t ; > f } 

i=0 i=0 i=0 i=0 

n—1 , n— 1 



(^ + 7,,=xp{E* + ('-I")e} 

1 1 H r < < f 'H f ' > f f e~~ 2d ^= s > ti dta...dt n -idt n 

^ 1=0 !=0 

(2^)"^" / exp{ £/^+ (f- X>W}l{ J> <Ad/ - 



<(2dye- 2d 'e^' n , _ (9) 

where the last line follows from Lemma 9. Using the definition of our function % y 
we get 



^'nA=^ few)+£), n 



< e m M ) + e)t e -n Y\ T( I T , , . (10) 

Next recall that p is fixed, and m n = \n p \ . Let 

G„ = {Z< 1 »,...,Z<;"»»}c{zeZ rf : | z |<n}, 

and call the complement G% the set of sites with very low potential. Note that there 
are at least | + \_8rt\ — m„ points in the path y that belong to G c n . Hence there are 
at least 

\Z^\+[Sr t \-m n 
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terms in the product in the left hand side of (10) that are smaller than 
provided this is less than 1. Combining this with (8), (9) and (10), we get 

u[ sM) (t) < £ e -'v( M «* ) + E - 2^/ )'(M^ ) -M^' ) )- |z "* )hL5 ''' J+m '' 

< (2d)"e-"e^ )+£ - 2 ^' (P^L\ogny^- ySrA+m " . 
Taking the log of the above and defining C £ := log(^) — log(£^) we get 

j log u[' M (t) < 2iogM +M W-2d + £-i(|Z»| + LSr,J -m n )log(^log«) 
<MW-i|ZW|loglog|Z«|-2rf + e + ?C e -5gg2+ (l), 

where we use that |Z„ | + [8r t \ < n. Observing that loglogn > (1 +o(l))loglogf 
and fC £ = o(l ), uniformly for all n in the given range, concludes the proof of (b). 

To prove part (a) we show that regardless of the distance travelled by the path, it 
hits a site with very low potential in every other step. Recall that a set H of vertices 
of If 1 is totally disconnected if there is no pair of vertices (x,y) € H 2 such that 
|*-y| = l. 

Lemma 10. Almost surely, for sufficiently large n, the set G„ is totally disconnected. 

Proof. We prove the statement for d > 2 first. If i and j are distinct integers in 
{1, . . . ,m„}, the random pair of points (Zn\z„ j) ) is uniformly distributed over all 
possible pairs of points in the ball of radius «. As no vertex has more than 2d neigh- 
bours, we have P(Z, l ! ' ) —Zip) < 2d/l n . Summing over all possible pairs i,j we get 

(m \ 2d 

P(G n not totally disconnected) < £p(Z« -Zf) < " — < Cn 2p - d . (11) 

tn 



KJ 



for some constant C. Since p < j and d > 2 we can apply the Borel-Cantelli lemma 
and obtain the result. We now prove the the same result when d = 1 . We introduce 
a new quantity 

m'„ = [n p '\ with p <p' <\ 

Let G' n be the set of the m' n vertices in the ball of radius « where the biggest values 
of £, are taken, and let p n be the biggest integer power of 2, which is less than «. 
Note that, by (1 1), the set G' Pn is totally disconnected for all sufficiently large «. 
We now prove that 

G n ^ = G' lpn for all sufficiently large n. (12) 
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For this it suffices to show that at least m„ points of G' 2pn are in the ball of radius n. 
Indeed, if we assume this and also that G„ <£. G' 2pn we can find a vertex zo satisfying, 
|zn| < n > Zo G G' 2pn and zo ^ G„. This implies that every zeG„ satisfies t, (z) > t; (zo), 
because G„ is the set where the largest values of ^ are achieved. Then, because 
zo € G' 2pn , we have 

G n c > §( Zo )}nB„ c {^( Z ) >§(zo)}ns 2p „ c g^„, 

which leads to a contradiction to our assumption. 

In fact we will prove the slightly stronger statement that there are at least ni2 P „ 
vertices of G' 2 in the ball of radius p n , and we will now write p instead of p n . We 
write 

G 2p = Wo,--- 

where d i is the vertex where £, (a-) = M^t 1 " and introduce 

* = (^)o<K^-i withX, = l {|a;| < p} and |X| = £ X*. 

Observing that m' 2p = o(p) and that G' 2p is uniformly distributed over all possible 
ordered sets and recalling that the box of radius p contains 2p+l vertices, it is easy 
to see that for p big enough, 

P(Xj = l\Xi=Xi,Vi <j)<\ m&P(X j = 0\X=x i ,'ii < j) < f, 

for all j < m' 2p — 1 and for all fixed (xq, ...,xj-i) G {0, 1} J . Hence 

P(\X\<m 2p ) = £ £P(X=*) 

= exp (-»4 p log(4/3) + (m 2p - l)\ogm' 2p + logm 2p ) 

= e -(2 P )P'(l +o( l))< e - n P' asn < 2/ , n . 

Using the Borel-Cantelli lemma we can prove (12), which implies the statement. 
We define the set of paths & n to be 

&n = fy = (yo, • • • ,y n ) ' yo = 0, \yi- 1 - yi\ = l , 

{y ,...,y„}n{Z^\...,Z^}^9}, 



so that 
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, n-1 



yE&>„ 



U["\t)= E (2rf)-»E[«p{£^(y,)+(t-£Ti)§(y„)} 

k i=0 i=0 

n-1 n . 



We can now argue similarly as for part (b) but using this time the fact that for any 
path in the number of step out of G„ is at least \n /2J . More precisely, 

U{"\t)< £ £-" e ( M ""+ £ - M >(M*' ) -M„ m " ) )- L ' ,/2J , 

and taking the log of the above and defining C £ := 21og(^) — log(£^) we get 

j log f/W (/) < 2 log M + M« - 2d + e - i L»/2J log ( ^ log n) 
= M<"-|r(loglog«-C £ )-2 £ / + e + (l), 
which concludes the proof of (a). 



2.4 Analysis of the variational problem 

We use the point process framework established in [6, Section 2.2] adapting the 
approach of [1 1, Chapter 3]. We only give an outline of the framework and sketched 
proofs here, see [6, Section 2.2] for more details. 

Observe that n{dy) := e~ y dy is a Radon measure on G := (— °°,°°\. For any 
z e If', ieM and r > 0, we have 

r d P{t, (z) - dlogr > x) = r d e- dlo s>- x = e~ x = /i ([*,«>]) • 

Define, for any q, x > the set H% := {(x,y) eR d xG: y> q\x\ + t}, where R d is 
the one-point compactification of W 1 . As in [6, Lemma 4.3] we see that the point 
process 

£-= E 5((f,§(z)-rflogr)) 

converges in law to the Poisson process £ with intensity Leb^ ® jU in the sense that, 
for any pairwise disjoint compact sets , . . . , K n C //? with Leb^ + i (<3/T) = 0, we 
have that (£ r (K\), £ r (K„)) converge in law to 

n 

0Poiss(Leb rf ® /*(£,■))• 

!=1 
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We further note that for z = f 1+0 W we have 

W,(z) := % (z) - ^loglog|z| = §(z) - M (i +o( i)). 

As in [6, Lemma 4.4] applied to T t (z,x) := (z, jc — |z| ) we infer from this the conver- 
gence of the point process 

°* : = L 5 ((^w(z)-^iogrO) 

in law to a Poisson process (U with intensity 

(Leb d ® At) o T," 1 = e-W-^dzdy, 

where now the compact sets K\ , . . . , K„ can be chosen from the set H z := W l+l \ 
(R d x (— oo, t)). The form of these and the previous domains, and in particular the 
use of the compactification, ensure that we can use these convergence results to 
analyse the right hand side of the final formula in Proposition 1 . 

Lemma 11. Let Xj- and X t be the sites corresponding to the largest and second 
largest value of \\f t (z), Z € Then ) — \p t QQ ) converges in law to a stan- 

dard exponential random variable. 

Proof. Using careful arguments in the convergence step we obtain, for any a > 0, 

P( Wt (Xl l) )- Vt (X t il} )>a) 

= x (y ,«)) =0,GJ,(r' x {y}) = l,®,(R d x (y-a,y))=0) 

y 

-> J P(®(R d x (y,oo)) =0)P(®(M. d x (y-a,y))=0)e~ y dy 
= j exp(-e- y+a )e- y dy = e~ a . 

Lemma 12. Let Xf be the site corresponding to the largest value of\\/ t (z), z € 
ThenX t w /r, converges in law to a random variable in R d with coordinates given by 
independent standard exponential variables with uniform random signs. 

Proof. As above we obtain, for any A c R d Borel with Leb c ((dA) = 0, 

P (*lL eA) =£p{a,{R d x (y.oo)) =0,(D t {A x {y}) = l) 
v r > 7 y 

-> J dzj dye-^- y P(®(R d x (y,oo)) = 0) 
= [ dz [ dy exp(-e~ y )e~ y ~^ = [ 2~ d e^dz- 

J A J J A 

Observe that this implies that the limit variable has the given distribution. 
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2.5 Proof of the almost sure asymptotics 



Note that combining Lemma 1 and Lemma 6 establishes the almost sure lower 
bound for the liminf result in Theorem 2. To find a matching upper bound, recall 
from Lemma 2 that, for sufficiently large t , 

U <N £ (t)-2d + e 

for N e (t) := Af c ( e )(f). We now approximate the distribution of N E (t). 

Lemma 13 (Approximation for the distribution of N s (t)). Let b t t °°> then 

log (P(N e (t) <bt))= -e- b >r*2 d (\+o{\)). 

Proof. Observe that 

P(N e (t)<b t )= n F^ + M(loglog| Z |-Q)). 

(/(Iog0 2 <|z|<'logf \ / 

The values which |z| can take are such that loglog |z| = loglogf + o(l) uniformly 
for all z, and since b, — > °o, we have, 

log (P(N e (t)<b t )) 

log(l-exp(-^-M(loglog?-C e + (l))) 

f/(logf) 2 <|z|<r logf 

= - (l + o(l)) £ e -b,-^+oW) 

r/(log() 2 <|j|<flogr 

= -e bt rf(\+o(\))j^e |A|(I+ " (I)) l{ loglog , /(log , ) 2< H < logfloglog ,} 



dx 



To obtain our final result, we apply the dominated convergence theorem to the inte- 
gral, which converges to 2 d . 

We are now ready to prove the upper bound. We consider a sequence of times 
t„ := exp(n 2 ) for which N e (t n ) are independent random variables, in order to use 
Borel-Cantelli. 

Lemma 14 (Upper bound for lower envelope of N £ (t„)). For any small c > 0, 
almost surely there are infinitely many n such that 

N e (t„) < d log t„ - ( 1 + d - c) log loglog/„ . 

Proof. Note that {N £ (t n )) n >N is a sequence of independent variables if is large 
enough. To see this it suffices to notice that the different (N £ (t n ))n>N depend on the 
values of the potential on disjoints areas. Indeed 
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> « 2 exp (« 2 ) = t„ \ogt n for all large n. 



t n +i _ exp(« 2 + 2« + l) ^ 



(log f , 7+1 ) 2 (» + l) 4 

Now we use Lemma 13 with b t = dlogt — (1 + d — c)logloglogf and we get, 

log(P(N e (tn) <b tn )) =-2 rf (loglogf„) 1 - c (l+ (l)) > -logn, 

for all sufficiently large n. Hence the sum over the probabilities diverges and we 
obtain our result by applying the converse of the Borel-Cantelli lemma. 



2.6 Proof of the weak asymptotics 

To prove Theorem 1 we show that the upper and lower bounds we found earlier for 
L t both satisfy the required limit statement. We first state the result of [6, Proposition 
4. 12], which describes the limit result for the lower bound N_(t). 

Lemma 15 (Weak asymptotics for N(t)). As t tends to infinity, 

N(t)-d\ogt +c/logloglogr ^>X, where P(X < x) = exp ( - 2 d e~ x ) . 

Next we check the analogous limit theorem for the upper bound N £ (t) and thus 
complete the proof of Theorem 1 . 

Lemma 16 (Weak asymptotics for N £ (?)). As t tends to infinity, 

N e (t) - dlogt + d\og\og\ogt =3>X, where P(X < x) = exp ( - 2 d e~ x ) . 
Proof. Fix x G R and apply Lemma 13 with b t = dlogt — d log log log? +x to get 

log (P(N £ {t) -dlogt +</logloglogf <x)) = -e- x 2 d (l+o(l)), 
which proves our result. 

2.7 Proof of the scaling limit theorem 

We recall thatx/* 1 (k = 1, 2) is the site at which 

= 7 i°gi°gkl 

takes its k th largest value. Fix 8 > and write 

U(t) = U 1 (t) + U 2 (t)+U 3 {t) + U 4 {t) + U 5 (t) + U 6 (t), 
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where U2 and U$ were defined in Lemma 7, resp. Lemma 8, and 

Ui (t ) = E [exp { jf 4 (X,.) d,} 1 { < J, < t a, , max 4 (X s ) > ' } 

l{X t W € {X s : < s < f}, |X f < <5r,} 

I/ 4 (0 = E [ exp ( f § (X s ) 1 { ta t < h < t log t , max | (X,) > M** 7 ' ' } 1 , 

L JO J 0<s<t 1 J 

f/ 5 (0 = E [exp { (X s )ds] 1 { < J, < t a, , max 4 (X s ) > 1 } 

<><*<*}}], 

t/ 6 (f) = E[exp{ jf$&)£fa}l{^ </, <fa„ max ><'''} 

l{X t w € {X s : < 5 < r}, |X, > 5r r } 



Observe that our result follows if the contributions of C/, (f ) for i = 2, ... ,6 to the 
total mass are negligible, as t/i(f) only contributes to the mass distributed on points 
close to X t { on the r t scale. 

Lemma 17. Suppose a t \, anc/ «r log log t — > 00. 77ie« we /1a ve, /« probability, 

Hm^=li m ^=lim^W=0. 

rt- £/(*) rt~ C/(0 rh- C/(r) 

Proof. For the first statement we use Proposition 1 (a) to see that 

-\ogU 4 {t)< sup \M^-^-(loglogn-C e )} + e-2d + o(l). 

1 n>ta, ^ M ) 

By Lemmas 12 and 11 the limit of the right hand side is strictly smaller than the 
growth rate of U(t), proving that the first limit in the statement equals zero. 

Using Proposition 1 (b) with 5 = and summing over all 1 < k < t a with ^ 
Zn \ and over all n with 2f/(logf) 2 <n<ta t we get 

-logt/ 5 (f)< max y t (z)-2d + o(l) = y/,(X r (2) ) - 2<f + o(l) in probability. 

By Lemma 1 1 we find e > such that, with a probability arbitrarily close to one 

j logf/s (t) < Wt (Xt m ) ~ 2d - e + o{ 1 ) , 

and a comparison with the lower bound N(t) for the growth rate of U(t) proves the 
second result. 
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For the third statement we use Proposition 1 (b) with the choice of 8 > from 
the statement. Summing over all 1 < k < t a and n with 2t / (log?) 2 < n < ta t we get, 
as above, 

jlogC/ 6 (0 < \lf t {X t m )-2d-S + o(l). 

We can now argue as before that this rate is strictly smaller than the lower bound 
N_(t ) for U (t ), proving the final statement. 

We can now complete the proof of Theorem 3. By definition we have 

1 > liminfv,{z 6 1 d I \z-X! l} \ < 8r t \ > liminf^-^ = 1 -limsupy ^$1. 

~ n» 1 1 1 " J " ft- 1/(0 , to e v p l U(t) 

Combining Lemmas 7, 8 and 17 we see that the limsup is zero, so that we get 
lim v t {z e Z d I \z - X, (l) I < Sr t \ = 1 in probability. 

(to" 

Combining this with the convergence of Xf' jr t given in Lemma 12 and recalling 
that 8 > was arbitrary concludes the proof. 



3 Concluding remarks 



It would be interesting to study scaling limit theorems for potentials with lighter 
tails and thus shed further light on the number of relevant islands in these cases. 

The techniques of the present paper appear suitable to treat cases where the rele- 
vant islands are single sites, which is the case for potentials heavier than the double- 
exponential distributions. For the double-exponential distribution itself and lighter 
tails, arguments related to classical order statistics of i.i.d. random variables need 
to be replaced by eigenvalue order statistics for the random Schrodinger operator 
A + E, on £ 2 (Z ), making the problem much more complex. Work in an advanced 
state of progress by Biskup and Konig [1] deals with the double-exponential case 
and strongly hints at localization in a single island of finite size in this and other 
cases of unbounded potentials. 

For bounded potentials the question of the number of relevant islands and the 
formulation of a scaling limit theorem at present seems wide open and constitutes 
an attractive research project. Sznitman in [12] discusses an 'elliptic version' of 
the Anderson problem, describing Brownian paths in a Poissonian potential condi- 
tioned to reach a remote location. Sznitman's technique of enlargement of obstacles, 
described in [14], offers a possible approach to the scaling limit theorem, leading 
in [13] to a study of fluctuations of the principal eigenvalues of the operator A + E, 
and moreover an analysis of variational problems somewhat similar to those that we 
expect to arise in the proof of a scaling limit theorem. 
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In the light of our result and this discussion it would be of particular interest 
to know whether there at all exist potentials which lead to more than one relevant 
island, and if so, to find the nature and location of the transition between phases of 
one and several islands. 
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